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1 Introduction 



In several publications [1-4] it is discussed that hadronic decay cannot be ignored in hadron 
models. As a consequence it is essential for the description of mesons (and baryons) to know the 
allowed decay modes and their relative strengths. 

The use of independent harmonic oscillators for the calculation of the angular-momentum part 
of rearrangement matrix elements, has been demonstrated in a previous paper [5]. There it was 
shown that the difficulties, which arise in naive recoupling schemes [6] when for a decay process 
the decay products have internal angular excitations, can easily be handled once one introduces 
explicit wave functions for the partons involved. This fact has already been elaborated in [7] 
and [8] and successfully applied to the calculation of branching ratios for the decay of mesons. 
Here we want to use this as a method for determining the coupling constants and transition 
potentials in a multichannel meson-meson scattering model [9,10]. 

There is, however, in principle no limit on the number of possible decay modes in the schemes 
of Refs. [7,8]. For this reason we prefer the rearrangement method of [5]. Essentially, when one 
wants to calculate all possible decay channels of a given initial state, then the harmonic-oscillator 
treatment has the advantage that there is only finite number of possible decay channels involved, 
which in a certain sense is complete. As a consequence, one can always easily verify whether all 
possibilities are taken into account. 

In this paper we will calculate explicitly the matrix elements for the decay processes: 

meson meson + meson 

(Sects. El El EI)- Of course, the Clebsch-Gordannery is straightforward, but it is useful to have 
the complete formulas written somewhere. Moreover, the here presented strategy is not totally 
trivial. We will see that in order to obtain the proper selection rules for the decay products, 
we only need to impose Fermi statistics on the quarks and the antiquarks. G parity is then 
automatically obeyed. In Sect. El some of the resulting coupling constants are compared to the 
data. The construction of a transition potential which can be used in models [9, 10] for the 
description of hadronic decay, is given in Sect. El 
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2 The meson just before decay 

Just before decay a meson is assumed to consist of four partons in the '^Po niodel (see also Fig.Q): 

1. the original quark- ant iquark (qq) pair which carries all the quantum numbers of the decay- 
ing meson, 

2. the newly created qq pair with the quantum numbers of the vacuum 



1 origvnalPf^^ 




Figure 1: The four-quark system of the original pair and the newly created ^Pq pair, just before decay. 

First we introduce a quarkspinor which describes a quark with flavorindex a (here the number 
of flavors is restricted to three; the generalization to more flavors is straightforward), colorindex 
a and spin magnetic quantum number n by 

qlif^) ■ (1) 

For an antiquark we use similarly the notation 

Qa^if^) ■ (2) 
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The wave function describing the system of four partons which compose a decaying meson, 
is the product of two parts: 



1. The wave function of the original qq pair, which consists of a color-singlet quark- ant iquark 
pair with (flavor, spin) indices (a, /ii) for the quark and {b, ^12) for the antiquark 

X^csoniJi Jz-J,s,n; a, b;r 12) = (3) 
E ' ^ Cj j ' 0n/,,,(ri.) -^q:{^^^)q:{^l2) ■ 

Here and in the following we adopt the standard summation convention for repeated indices. 
The antiquark is at relative position ri2 with respect to the quark. This system has relative 
angular momentum £, total spin s, total angular momentum J {z component: Jz) and radial 
excitation of the spatial relative motion n. 

2. Similarly the newly created color-singlet and fiavor-singlet ^Pq pair is described by 

^3Po(r34) = E \ ° c j j \ </'o,i,M(r'34) ^ g^^a) (/i4) • (4) 

where we have taken the lowest radial excitation in the spatial part. The precise definition 
of the harmonic oscillator wave functions <pn/.m (r) can be found in [5]. 

The total wave function for the four particles must be antisymmetric with respect to the 
interchange of either two quarks or two antiquarks. For this purpose we define the exchange 
operator which operator interchanges partons i and j. This way we obtain for the wave 
function for a meson just before decay 

;i_pi3_p24^pi3p24^ \M + ^Po; J,Jz;i,s,n;a,b;ru,ru,ri2,34) , (5) 

where 

\M + ^Po; J, Jz] i, s, n; a, b; ri2, r^, 12,34) = 

= 'Xmcson{J,Jz-J,s,n;a,b;ri2)X:ip^{ru)(j)o,o,o{ri2,34) ■ (6) 

In Eq. ©, the relative motion of the two qq pairs, which have relative position ^12, 34, is assumed 
to have the ground state quantum numbers. So far we have not taken into account normalization 
factors, but we will come back to this point in Sect. |31 
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3 The meson+meson final state 



When, out of the two qq pairs, two new mesons are formed, we have a system of two distinct 
objects, each of which consisting of a color singlet qq pair (Fig.|2|). We will not treat here possible 



q {d', ^i'^) 




Figure 2: The four-quark system of the original pair and the newly created ^Pq pair, just before decay. 

octet qq pairs. The quantum numbers (total angular momentum, orbital angular momentum, 
spin and radial orbital excitation) of the mesons are given by ( J^, £a, Sa, ria) for meson 1 and by 
{Jf,,ib,Sb, rih) for meson 2, whereas their flavor indices are given by (a', d') and (c', b') respectively. 

The wave function must be symmetric under the interchange of the two mesons, or in terms 
of the quarks, symmetric under the simultaneous interchange of quarks and antiquarks, i.e. 



^_^pi3p24j |7\^_^jvf . jj^.^ Jc,ic,nc; Ja,L,Sa,na; Jb,ib,Sb,nh; a' , d' , c' ,b'; r-^, rga, 7-14,32) 



(7) 



where 



\M + M; J,J^; Jc,ic,nc; Ja,ia, Sa,na; Jb, h, Sb,nb; a',d',c',b'; r^, r32, ri4,32) = 

^ Jc J /^Ja Jb Jc y^^a Ja ^^b ^b Jb 

~ ^ M,m,J, MaMbM, rriafiaMa rUbfibMb 

ma,fJ.a 

nil,, fJ-b 

(Pna,la,ma (^u) (pn,A,m, {t 32) 0n./.,mc (r" 14,32) E \ la (A^'i) Qd' if^-' a) Qp (Z^'s) Qb' if^'' 2) 

I I "J 

M 1:M 2 
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4 The matrix elements 



We obtained now two wave functions of the following structure: 

1. before decay (see Eq. Eland Fig. IH): 

2. after decay (see Eq. [7|and Fig. El): 

(l + P^^P^^) \M + M) , (10) 

We demand the wave function expressed by the ket |M + M) in Eq. (jiup to be such that the 
following matrix elements vanish, 

(M + 3Po|pi3|M + M) = {M + ^Pq\P'^^\M + M) = , (11) 

i.e., when we interchange only the quarks or only the antiquarks, then the wave function must 
have such a structure that the result is orthogonal to the states we consider. This has to do with 
the fact that nature seems to know which quark and antiquark belong to each other. As we may 
see, it also expresses the fact that nature obeys the OZI rule, or satisfies a kind of string idea. 
We ignore here the fact that, of course, in nature the relation (fTT|) will only be approximately 
satisfied, because the meson wave functions have some spreading. The here allowed processes are 
(OZI allowed) 

ab —>■ aq + qb or ab qb + aq , (12) 
the non-allowed processes are (OZI forbidden) 

ab —>■ ab + qq or ab ^ qq + ab . (13) 

The matrix-element 

M = (M + M|(l+Pl3p24)t^l_pl3_p24^pl3p24)|^^3p^^ ^ (14) 

involves the product of exchange operators 

'l + Pl3p24y _ pl3 _ p24 ^ pl3p24\ ^ 3 (l - P^^ _ p24 ^ pl3p24\ _ (^5) 



When we combine Eq. (fT^l) with the conditions (fTTj) . then, besides a factor 2, we obtain for the 
transition matrix element (fT^ : 

M = {M + M;rii,rs2,rM,32\{l + P''P^')\M + ''Po;ru,ru,ruM) 

= {M + M;ru,r,,,ru,34\iP'' + P^^)\M + ^Po;ru,r^s4,ru,34) ■ (16) 



Since the normalization constants are not yet relevant, we prefer to take formula (jl6j) for the 
definitions of the wave functions, rather than the wave functions defined in Eqs. (0) and ((Tj), i.e. 



1. before decay: 



2. after decay: 



M [-/, J.; s, n- (a, b)] [P'^ + |M + , 



(17) 



M2[J,Jz\ Jc,ic,nc; Ja,ia,Sa,na; Jb, ib, Sb, rib, {a',d'),{c',b')] \M + M) . (18) 

When we use the harmonic oscillator wave functions which are defined in [5], then we obtain 
for the normalization constants in Eqs. (jl7|) and (fTBj) . 



J\fi [J, Jz] i, s, n; (a, b)] 



+ Sjfi6j^^oSe,i6s,i6n,o (ab\{l} 



(19) 



which is only different from 1 when the original quark-antiquark pair (a, b) is in a flavor- singlet 
^Pq ground state. P^^ interchanges the quarks and the antiquarks. The bracket {^afe|{l}^ 
represents the coefficient of the SU (3)-singlet component in the full flavor wave function of the 
initial meson. 

For N'2 we find 



M2[J,Jz] JcJcUc] Ja,L,Sa,na; Jb, ib, Sb, rib] {tt' , d') , {c , b')] 



(20) 



The spatial parts of the matrix element |T6|) consist out of the diagrams as represented in 
Fig. |31 (see Ref. [5] for the definitions). The internal lines of diagram Q are function of the 
matrix elements an, a^^ of a 3 x 3 matrix a. When we restrict ourselves to the case of equal 
quark masses, then a is for the first term of (jl6|) given by 

r \ 



/ ^32 \ 

V ^'32,14 / 

and for the second term by 

/ \ 

V »^14,32 J 



I 1 

2 

1 
2 



/ 1 

2 

1 
2 



V 2 



/ ^12 \ 

V ^12,34 ) 

( ri2 \ 

V ^^12,34 J 



(21) 



(22) 



7 



{n,£, III) 



(0,1, m) 



(0,0,0) 



52 - 



_X OL22 —> 



^ ^^^^ ^ 



9- 

^33 — -"^ 



Figure 3: Diagrammatic representation of the spatial parts of formula p6() . 

We denote the respective results of the diagram of Fig. El for the a's defined in Eqs. (j2l|l and 
respectively by 



di {n, £, Hi] fi; Ua, 4, ma, m, 4, mb, Uc, 4, rric) , (23) 

and 

d2 {n, i, jj.; Ua, 4, ma, Ub, ib, mb, n^, 4, m^) . (24) 
With these definitions we obtain for the matrix elements pfij) : 

M = A/-, I J. n; (a. 6)] ^ { (25) 

"la , Ma , , /ij, 

mc,/if, /is, Ml 1/^2 
mi , m2 , ms , m4 

niaHaMa uib ^b Mb jJ-i jJ-s Jz AH At2 mi 7712 m^m^Hb 

i 1 1 S i i 1 

<^ di {n, i, Hi, /i; ria, 4, "^a; ^fe, 4, mb, n^, ic, C J J C J J ,5^,^, ,5^,^ + 

7773 7772 /^s '^l 

i i S i i 1 1 

+ 4 (77,£, /ii; /i; na,ia,ma; nb,ib,mb; nc,ic,mc) C*^^^^^ ^11137112^2 ^'^''^ '^'^'^"^'''M ' 

In Sect.Elwe smuggled two assumptions into the procedure, which had to do with the excitations 
of the ^Po-pair. This resulted in the four zero's in diagram (jS)). The generalization to other 
quantum numbers is straightforward but will not be treated here. 
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5 Results and comparison with experiment 



In this section we study the results of expression for some quark- ant iquark systems. First 
we introduce the nomenclature which has been used here for mesons. In Table ^ we give the 
precise meaning of the symbols which we use in Table El 



meson 


qq 




jPC 


poc U-U-UoL^didi D 






n~+ 

U 




first radial excitation 


2' So 




vectors 


ground state 


I'Si 


1— 




first radial excitation 


2' Si 






second angular excitation 




I"— 


scalars 


ground state 


I'Po 


0++ 


axial vectors 


ground state {S = 1) 




1++ 




ground state {S = 0) 


I' Pi 


1+- 



Table 1: Nomenclature of quark-antiquark systems. The columns contain respectively, the most current 
denomination for mesons, the assumed qq state, the g^g-quantum numbers (n is the assumed radial 
quantum- number, S represents the total spin, L the orbital and J total angular momentum), and the 
more common quantum numbers {P stands for parity and C for charge conjugation; excitations are 
indicated by accents). 
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In Table 121 all [i.e. within the above specified assumptions) possible decay channels for the 
lowest radial and angular excitations of pseudoscalar, vector and scalar mesons are given. Note 





0-+ 


1— 


0++ 


(0-+ 


,0-+) 




-± 

24 






(0-+ 


,0'~+) 










(0"+ 


,1"-) 






(l^P) 






,0++) 










(0-+ 






1 

12 




1 (1^^) 


(0-+ 






1 

24 






(1- 


,1"-) 




-± 

72 


(iip) 






,1") 






(i^p) 














Ti4(l'^) 


(^ — 


1 // — \ 

1 J- ) 










(1- 






1 
6 










1 /-I 1 


1 

12 




i P) 


(0++ 


,0++) 










(0++ 


,1-) 




1 
8 






(1++ 












(1^" 


,1^-) 











Table 2: Three- meson vertices |(MM|y|M)p. The signs in front of the matrixelements are explained 
in the text. The orbital quantum numbers of the meson-meson system are given between brackets, in 
the spectroscopic notation v?^~^^L, where n is radial quantum number, S the spin and L the angular 
momentum. 

that the figures in each column of Table El add up to 1. This expresses the fact that each wave 
function (fT7|) can be completely decomposed into the basis (fTH|) . because both sets are solutions 
of the same Hamiltonian. 

The signs which are given in Table |21 are the results of the ratios of the terms stemming from 
P^^ and the terms stemming from P^^ in formula ()l(i|l . These signs are essential for conservation 
of G-parity, which can easily been seen once the 5'?7(3)-fiavor content is taken into account. 

We are now prepared to test the results of Table |21 to the available data. 
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In general the width of a decaying particle is given (see e.g. Ref. [11]) by 

T(m2, cos(i9)) 



Arrr 



dcos{'&) 



(26) 



The matrix element T in Eq. (j^Bj) includes the calculation of the overlap between the free meson 
wave function and the harmonic oscillator basis, which we have chosen. In fact we should do 

{MM\V\M) = {MM\n') {n'\V\n) {n\M) , (27) 

n,n' 

where \M) is a specific initial state of the qq mesonic system and \MM) represents the final state 
of the two decay products. But, since decay widths can be better calculated from the solutions of 
the multichannel Schrodinger equation involving the transition potential (this is done by us in a 
different program [9,10]), we will not be too rigorous in demonstrating that the matrix elements 
()26p make sense. So we assume that {n\M) is only nonzero for one value of n, and that then 
also the corresponding {MM \ n') is the only nonzero term. The overlap of the plane waves for 
the two meson system and a harmonic oscillator contains a momentum dependence of the form 

T (X . (28) 



where i is the relative orbital angular momentum of the decay products. Relation ()28|) . together 
with the phase space factor, gives 

r oc . (29) 

So, we arrive at 

r oc k^'+^ \{n'\V\n)f , (30) 

where n represents the radial quantum number of the decaying meson. A more precise form of 
formula (jHUj) is given by the imaginary part of the second order perturbation contribution to the 
energy within the model [9,10], as is given in [12,13]. We will, however, take as sufficiently 
accurate for our purpose here. 
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As examples we study the decay widths of e(1300) (nowadays /o(1370)), a2(1320), and 
Kq{14:30), which particles have enough two particle decay modes, to make a comparison of 
relation to the available data possible. The resuhs are given in Table El 



processes 


Eq. dnni) 


PDG80 


PDG04 


KK 

f ^ TTTT 


0.22 


0.11 


0.12±0.06 - 0.91±0.20 


a2 — > PTT 


(1 37 


20 - 25 


213+0 020 


K* Ktt 
K* K*n 


2.28 


1.6 - 2.1 


2.02±0.14 


K* Kn 
K; -> Kp 


7.76 


4.3 - 6.1 


6.7±1.1 


K* Kr^ 
K* Ku 


2.93 


- 3.7 





Table 3: Some branching ratios of e(1300) (nowadays /o(1370)), a2(1320), and K^{U30). Results of 
expression H3U|) compared to experiment: PDG80, C. Bricman et al., Rev. Mod. Phys. 52, 1 (1980), and 
PDG04, S. Eidelman et al., Phys. Lett. B592, 1 (2004) 



12 



6 The transition potential 



Inspired by the reasonable results of the comparison of our handwaving calculations with the data, 
we will derive in this section a transition potential for a coupled channel Schrodinger equation 
as described in [9,10]. In [5] we have shown how we might arrive at a local approximation of the 
potential for transitions from a quark-antiquark channel to a meson-meson channel. The outline 
given there will be made more precise in this section. 
The potential is given by 



where {i, ne) are the orbital quantum numbers of the qq mesonic system, where m is the quark- 
mass and where u is the universal oscillator frequency, and by 



where {£', n'^) are the orbital quantum numbers of the two meson system. 

The matrix elements {n' \ V\ n) are those given in [5]. For the lowest lying qq systems n' is 
linearly related to n. 



{MM,x\V\M,x) = {MM,x\n') {n'\V\n) {n\M,x) 



(31) 



The wave functions are in lowest order assumed to be given by 




(32) 




(33) 
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Figure 4: Transition potentials for the decay of pseudoscalar mesons (r and V{r) are dimensionless) . 
The curves are determined with expression H31|) for the five three- meson vertices of Table @ . 
The upper curve (blue) corresponds to the pseudoscalar vertices with (vector, vector) and (pseudoscalar, 
vector) the middle curve (violet) to the pseudoscalar vertex with (vector, positive-C-parity axial vector), 
the lower curve (red) to the pseudoscalar vertices with (pseudoscalar, scalar) and (vector, negative-C- 
parity axial vector). 

In Fig. m we present the results of Eq. (j31|) for the decay processes of pseudoscalar particles. 
Although in Table El five different types of decay channels for a pseudoscalar meson are shown, 
we find that the spatial behavior for two pairs is exactly the same. We also observe from Fig. |3] 
that the decay processes P — > PV and P — VV are dominant. This result is important because 
it justifies the procedure of [9, 10]. There we have limited us to these channels. For doing so we 
now have two arguments: the thresholds are low, the couplings are largest. We also find that the 
form of the here calculated transition potential is in agreement with the one from [10]. However, 
in [10] the position of the peak is a parameter which can be fitted to the data, whereas here it is 
a result of our procedure. 
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Figure 5: Transition potentials for the decay of vector mesons (r and V{r) are dimensionless) . The 
curves are determined with expression for the nine three- meson vertices of Table (jSJ. 
The curves, from the curve with the highest maximum to the curve with the lowest maximum, represent 
the vector vertices with respectively (vector, vector in P wave, total spin 5 = 2), (pseudoscalar, vector), 
(vector, positive-C-parity axial vector), (scalar, vector), (pseudoscalar, pseudoscalar), (pseudoscalar, 
positive-C-parity axial vector) and (vector, negative-C-parity axial vector), (pseudoscalar, negative-C- 
parity axial vector), (vector, vector in P wave, total spin 5 = 0). 

In Figs. El and ini are depicted similar results for the vector and the scalar mesons respectively. 
Is it reasonable in the case of pseudoscalar mesons to select two dominant decay modes, in the 
case of vector mesons the situation is less comfortable (see Fig. Ej) and becomes even hopeless in 
the case of scalar mesons (Fig. IHl). 

So, in the next stage of the program of the multi-channel model of [9, 10] all possible decay 
channels should be taken into account in order to decide afterwards which of them are important 
for the description of a specific hadron. The here presented method can supply us with the precise 
form of the transition potentials. Finally we must notice that the potentials of the Figs. |3J Eland 
El exhibit a peaked structure. In [12] it is argued that such types of potentials can very easily 
explain the radial hadron spectra, especially the large level splittings between the ground states 
and the first radial excitations with respect to the other mass differences. 
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I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

0.5 1.0 1.5 2.0 2.5 , 



Figure 6: Transition potentials for the decay of scalar mesons (r and V{r) are dimensionless) . The 
curves are determined with expression (|31j) for the eleven three-meson vertices of Table ((2j. 
The curves, from the curve with the highest maximum to the curve with the lowest maximum, represent 
the scalar vertices with respectively (vector, vector in D wave, total spin 5 = 2), (pseudoscalar, positive- 
C-parity axial vector) and (vector, negative-C-parity axial vector), (pseudoscalar, pseudoscalar), (vector, 
vector (angular excitation)), (pseudoscalar, pseudoscalar (radial excitation)) and (negative-C-parity 
axial vector, negative-C-parity axial vector), (scalar, scalar), (vector, vector in S wave, total spin 5 = 0), 
(vector, vector (radial excitation)), (positive-C-parity axial vector, positive-C-parity axial vector). 
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Conclusions 



We have shown here that the method which was roughly outhned in [5] , can indeed be used for the 
construction of transition potentials to be used in coupled channel models [9, 10] for mesons (and 
baryons). The procedure is still not unambiguous, but, for a first inspection of the results, very 
suitable. The next refinement could be, to take different quark masses in the matrix elements 
This leads to different transition potentials for decay due to the creation of up and down 
quark pairs as for decay due to the creation of strange quark pairs. Another possibility is, to 
repeat the whole procedure for the 5*0(3,2) wave functions of [14], which are probably more 
realistic than harmonic oscillator wave functions. 
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